Non-Gaussian Density Fluctuations from Entropically Generated Curvature 

Perturbations in Ekpyrotic Models 
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We analyze the non-gaussian density perturbations generated in ekpyrotic/cyclic models based 
on heterotic M-theory. In this picture, two scalar fields produce nearly scale-invariant entropic per- 
turbations during an ekpyrotic phase that are converted into curvature modes after the ekpyrotic 
phase is complete and just before the big bang. Both intrinsic non-linearity in the entropy pertur- 
bation and the conversion process contribute to non-gaussianity. The range of the non-gaussianity 
parameter fNL depends on the details of the scalar field potential during the ekpyrotic phase, and 
on how gradual the conversion process is. Although a wider range is possible, in principle, natu- 
ral values of the parameters of the potential combined with a gradual conversion process lead to 
values of —60 J; /jvl +80, typically much greater than slow-roll infiation but within the current 
observational bounds. 



Ekpyrotic [1] and cyclic [2] models of the universe use 
quantum fluctuations of scalar fields produced during a 
slowly contracting phase with equation of state w > 1 
to generate the observed nearly scale-invariant spectrum 
of curvature (energy density) fluctuations after the big 
bang. One mechanism considered for converting the fluc- 
tuations of a scalar field into cosmological curvature per- 
turbations relies on higher-dimensional effects at the col- 
lision between orbifold planes (branes) along an extra di- 
mension [3]. Recently, however, a new "entropic" mech- 
anism [4, 5] has been proposed that relies on two (or 
more) scalar fields and ordinary 4d physics, stimulating 
new approaches to ekpyrotic and cyclic cosmology that 
may not require branes or extra dimensions at all [6-9] . 

In this paper, we wish to consider an important 
byproduct of the entropic mechanism: a non-gaussian 
contribution to the density fluctuation spectrum that is 
more than an order of magnitude greater than for con- 
ventional inflationary models and that can satisfy current 
observational bounds, including the recently claimed de- 
tection of non-gaussianity [10]. Our results differ signifi- 
cantly from the cases considered by Koyama et a/. [11] and 
Buchbinder et aZ. [12] in which they assumed an ekpyrotic 
(ui ^ 1) phase that continues all the way to the conver- 
sion of entropic to curvature fluctuations, as in the "new 
ekpyrotic" model [6]; for these cases, the non-gaussianity 
is amplified by the non-linear evolution on super-horizon 
scales to the point where Jnl, the parameter character- 
izing the non-linear curvature perturbation [13], reaches 
magnitude 0{c{), where ci « 2 -^/uT+T parameterizes the 
steepness of the scalar field potential during the ekpyrotic 
phase. A potential problem is that a minimum of ci > 10 
is required just to satisfy the current upper bound con- 
straints on ris , the spectral tilt of the scalar (energy den- 
sity) perturbation spectrum [5] ; and ci > 30 is needed to 
reach the best-fit value Ug « 0.97. Yet, excluding finely- 
tuned cancellations, the resulting value of /nl obtained 
in Ref. [6-9] is marginally consistent with current obser- 
vational bounds on non-gaussianity only for ci restricted 
to a narrow range (< 15). 

Here we show that, by having the conversion occur af- 



ter the ekpyrotic phase, as is natural from the point of 
view of heterotic M-theory and the cyclic model [2, 5], 
the non-gaussianity is reduced and much less sensitive to 
ci, so that one can naturally fit the best-fit value of 
and current bounds on both /jvl at the same time. More 
specifically, the net value of /tv l depends on both the in- 
trinsic non-linearity in the entropy perturbation as it is 
generated during the ekpyrotic phase as well as the dura- 
tion of the conversion process. These two contributions 
to f^L compete with each other, in that the intrinsic non- 
linearity, which depends on the steepness of the scalar 
field potential, is always positive and the contribution de- 
pending on the duration of the conversion process is neg- 
ative. In principle, with arbitrarily steep potentials and 
arbitrarily sharp conversions, a wide range of /jvl is pos- 
sible, and one can even have a finely-tuned cancellations 
of two large but oppositely signed contributions that give 
nearly zero non-gaussianity. Excluding these fine-tuned 
cases, though, we find that /nl lies within the current 
observationally favored range of 26.9 < Jnl < 146.7 that 
was recently reported at 2a [10] or the bound obtained 
earlier by the WMAP coUaboration, -36 < Jnl < 100 
[14], while, at the same time, matching the best-fit value 
of Us allowed by WMAP. 

As a concrete example, we will analyze the 4d effec- 
tive field theory derived from the ekpyrotic/cyclic model 
of [5]. The 4d theory describes the phase when two 
boundary branes approach one another and collide [15] 
in heterotic M-theory [16]; for a short review see [17]. 
The brane worldvolume is described by gravity and two 
canonically normalized scalar fields, (/>i,2, which param- 
eterize the distance between the two branes as well as 
the volume of the internal Calabi-Yau space [15]. For 
the purposes of this paper, though, the connection to M 
theory is not essential; it simply provides physical mo- 
tivation for considering this particular 4d field theoretic 
example. The same effects will occur in other 4d field 
theories in which the conversion of entropic to curvature 
perturbations occurs after the ekpyrotic {w 3> 1) phase. 

We assume a flat Friedmann- Robertson- Walker (FRW) 
background with line element ds"^ — —St^ + a?{t)5x^ and 



2 











oa 






^after / 










^-s,,^ before 




/ 




8s ^ \ 





fact, one could even imagine that it is the effective po- 
tential associated with this reflection which causes the 
ekpyrotic phase to end. In any case, the conversion is 
followed by a period of scalar field kinetic energy dom- 
ination, which sets simple cosmological conditions near 
the big crunch/big bang transition [2]. 

It is useful to recast the evolution in terms of the adi- 
abatic variable a (pointing along the background trajec- 
tory) and the entropy variable s, pointing transverse to 
the trajectory (see Fig. 1) [18]. Up to unimportant addi- 
tive constants that will be fixed shortly, we define 



FIG. 1: The background trajectory in scalar field space is 
straight everywhere except for a brief period after the ekpy- 
rotic phase around the time of reflection. The shaded region 
indicates the geometrically forbidden region bounded by the 
ajcis (/)2 = 0. In the presence of brane-bound matter, the tra- 
jectory reflects smoothly off this boundary. Also shown are 
the directions of adiabatic (5(t) and entropic (is) fluctuations, 
as well as the angle of the trajectory before and after the re- 
flection. 



scale factor a(i), and steep, negative, scalar field poten- 
tials of the form 



1^(01,02) = -1/16-^^1^^1-^26-^"^^ 



(1) 



where ci = ci(0i), ci = 02(^2) and Vi,V2 are positive 
constants. During the ekpyrotic phase, the ci are nearly 
constant and Ci ^ 1 . Then the Einstein-scalar equations 
admit a scaling solution 



-ln(-Vc^V./2t), 



1 _ ^ 2 



(2) 

which describes a very slowly contracting universe with 
e ^ 1 and w + X — 2t/'i':^ 1- the defining characteristics 
of an ekpyrotic phase. At the end of the ekpyrotic phase, 
the effective values of Ci ^ and the ekpyrotic potential 
energy ceases to be important cosmologically. 

The colliding branes solution corresponds to an almost 
everywhere straight line trajectory in scalar field space 
in the 4d effective theory. The only deviation from this 
straight trajectory is one required by the fact that the 
scalar field moduli space admits a boundary at 02 = 0, 
which from the higher-dimensional point of view corre- 
sponds to the locus where both the scale factor on the 
negative-tension brane and the volume of the Calabi-Yau 
manifold reach zero. As shown in [15], the trajectory 
simply gets reflected at this boundary. In the presence 
of matter on the branes, the reflection is smoothed out 
by an effective repulsive potential [15] and the trajec- 
tory never actually touches the boundary, as shown in 
Fig. 1. This smooth reflection is relevant for the den- 
sity fluctuation spectrum, because it suffices to induce 
the conversion from entropic to curvature perturbations 
without introducing any extra elements to the theory. In 
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where a = \J cpi + {cr \s positive by definition) . A dot 
represents a derivative with respect to FRW time. Then 
we can expand the potential up to third order as follows 
[20]: 



Vet 
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[l + -^ + ^^^/V], (5) 



where K3 is of 0(1) for typical potentials (the case of 
exact exponentials corresponds to K3 — — 4a/2/3). The 
scaling solution (2) can be rewritten as 

a{t) = {-tY/^ o- = -y^ln (-v/eV^i) s = 0. 

(6) 

The perturbations about this contracting universe have 
a growing mode - the entropy perturbation ~ corre- 
sponding to the relative fluctuation of the two fields: 
5s = 54>2 — 4>2 5<pi)/ a, We will decompose the entropy 
perturbation into a linear, gaussian part and a second- 
order perturbation by writing Ss — 6s^^^ -|-(5s^^^. Its equa- 
tion of motion, on large scales and up to second order in 
field perturbations is then given by [18] 

,5s + 3HSs +(Vss+ SeA 6s 



a 



a 2 



la a ] a 



Here H = a/a is the Hubble parameter, Va denotes a 
derivative of the potential along the background trajec- 
tory and Vs---s denote successive derivatives along the 
5s direction. The angle 9 of the background trajectory 
is defined as in [19] by cos{9) = ^i/ct, sin(6') — (f)2/cr. 
A useful expression for the time variation of this angle 
is ^ = —Vs/cr. In the simplest models, the background 
scalar field trajectory is well-approximated as a straight 
line [15], = 0. In this case, we have 02 = 70i, which 
implies ci = 7C2. For the colliding branes solution with 
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empty brancs, 7 = -^OH{t — trej)- where 9h is the Heav- 
iside step function (with values ±1 for t ^ iref) and tref 
denotes the time of the reflection, see [15] for details. The 
last term in equation (7) is a non-local term proportional 
to the difference in spatial gradients between the linear 
entropy perturbation and its time derivative. This diflfer- 
ence evolves as a^^ [18], so that it remains approximately 
constant during the ekpyrotic phase when a is very slowly 
varying. This ends up being exponentially suppressed 
compared to the entropy perturbation itself, which grows 
by a factor of 10^*^ or more during this same period [5]. 
After the ekpyrotic phase has ended, the non-local term 
grows because a oc {—t)~^/^, but this growth is negligi- 
ble compared to the exponential suppression during the 
ekpyrotic phase. Hence, the non-local term can be safely 
neglected. For analyzing the equation of motion, it is con- 
venient to use conformal time r. Denoting r derivatives 
with primes, and introducing the re-scaled entropy field 
5S = a(T) 5s, to first order in perturbations, equation (7) 

reduces to 5Sl" + (^^ - ^ + 5Sl = 0. The fast- 

roll parameter e = &'^/{2H^) « c^/(2(l + 7^)) defines the 
equation of state of the background scaling solution in the 
ekpyrotic contraction phase. In terms of e and its deriva- 
tive with respect to the number A/" = \n.{aendHendl o,H) of 
e-folds remaining before the end of the ekpyrotic phase, 
the solution to the equation of motion at long wave- 
lengths and linear order is 

5SL = fL{k,e){-T)-^+^/'-^. (8) 

Since e oc {w + 1) ^ 1 during the ekpyrotic contrac- 
tion phase, this leads to a nearly scale-invariant spec- 
trum with index [5] = 1 -I- | — ^j^- As discussed in 
[5], achieving a red spectral tilt consistent with best-fit 
measurement (rig « 0.97) typically requires e ~ J\f^-^~^ 
or ci > 30. 

On large scales and at linear order, the comoving cur- 
vature perturbation evolves according to [19] TZ = ^9Ss. 
Thus, the entropy perturbation generates a curvature 
perturbation when the background trajectory bends at 
the reflection. Moreover, this curvature perturbation is 
of the right amplitude if the potentials turn off at a time 
tend ~ — 10^ Mpj^, where Mpi is the Planck mass [5]. 

Next we consider the non-linearities responsible for the 
non-gaussianity of the perturbation spec;truni. The in- 
trinsic non-gaussianity in the entropy perturbation is pro- 
duced during the ekpyrotic phase and can be determined 
from the equation of motion for the entropy field to lead- 
ing order 

5S" +(^k^-^^5S+ ^aVsss {SSf = 0. (9) 
The last term is approximately 

The solution to the equation of motion (9), at long wave- 
lengths and up to second order in field perturbations, is 



then given by 



SS = SSL + ciSSLf, (11) 




So c characterizes the intrinsic non-linearity in the en- 
tropy perturbation. As a reference, we note that for 
empty branes, and with the original potential (1), i.e. 
with constant ci = — C2/-\/3, we have K3 = —4-^2/3 
hence c = — -\/e/6 = — ci/4. 

The time evolution of the curvature perturbation to 
second order in field perturbations and at long wave- 
lengths is given in FRW time by [18]: 

n = ^96s +49^){5s^''^f + ^6s6s]. (13) 
a a a 

At linear order a bending (6 ^ 0) of the background tra- 
jectory combined with a non-zero entropy perturbation 
source the curvature perturbation on large scales and re- 
sult in a linear, gaussian curvature perturbation 

Tli^^ J^eSs^^). (14) 

During this process of conversion, the entropy perturba- 
tion evolves according to Eq. (7) in a rather complicated 
way, which is why the above integral and the following 
calculation have to be performed numerically. Qualita- 
tively, a sharp reflection leads to a drastic diminution 
in the entropy perturbation and a much smaller value 
of the above integral, while a smooth and more gradual 
reflection results in an enhanced efficiency of conversion. 

The leading non-gaussianity is generated when modes 
are outside the horizon, so the non-gaussianity is of the 
local, wavelength-independent type [13, 21, 22]. The non- 
linearity parameter /jvl is then defined by the relations 
[13] = + /nl^l^ where is Bardeen's space- 
space metric perturbation [23] and is its value to lin- 
ear order. We adopt the convention of Ref. [18] where 
the comoving curvature perturbation is equal but oppo- 
site in sign to Bardeen's ( variable on large scales, so 
that during matter domination we have TZ = — |$jif , and 
thus TZ = TZl — ^fNhT^L- (Note that this definition of 
Inl agrees with Komatsu and Spcrgel [13] and WMAP 
[14] but is opposite in sign to Maldacena [22].) 

At the quadratic level, Eq. (13) then implies that there 
are three distinct ways in which the curvature perturba- 
tion can acquire non-linear contributions: 

(1) Before the reflection, the terms proportional to Vss 
and V^a create non-linearities in TZ during the entire ekpy- 
rotic phase in which the entropy perturbation is gener- 
ated; during this phase, the background trajectory is a 
straight hue, so Eq. (13) reduces toTZ= ^(— V"ss((5s)^ -|- 

^(5si5s), which to leading order reads TZ = ^/K— i)~^- 
The integrated contribution thus amounts to 

TZintegrated = ~^{^^ind)'^ ' i^^) 
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This term leads to 
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/integrated 
NL — 



(16) 



and is typically small, of order 0(10) or less. 

(2) At the reflection, any intrinsic non-linearity in 6s 
will be converted into a non-linearity in TZ due to the first 
term in (13); thus we get 



^intrinsic 

JNL — 



37^i 



(17) 



This term is positive and increases with increasing c, al- 
though generally not quite as fast as linearly because of 
the non-linear evolution of Ss'-^\ 

(3) Also at the reflection, the terms involving (^s*-^-* lead 
to non-linearities in TZ according to 



pre flection 

JNL ~ 



37^i 



-4^2)(5s(l))2 



ref 



a 



(18) 



In order to close the equations for the purposes of this 
study, we are modeling the reflection in terms of an effec- 
tive potential V that acts when 02 gets small and causes 
it to reflect; we find that the results do not depend sen- 
sitively on the form of V except for how gradual the re- 
flection is. Even disregarding the VJj^ prefactor, this 
term generically increases in magnitude as the reflection 
is made sharper and it is always negative in sign. 

The total /jvl is the sum of all the above contributions. 
In Fig. 2, we have shown the results for a typical value 
e = 36 (corresponding to e.g. c\ = C2 = 12) and for the 
range — 5 ^ K3 ^ 5 with the lower curves corresponding 
to K3 = — 5 and the higher curves to K3 = -1-5. In order to 
see how robust the results are, we experiment with dif- 
ferent potential forms to model the reflection of 02- The 
results shown here are for the potential forms Vi((/)2) — 
?;(02^+?''/'^^) andF2((/)2) = -y[(sinh(/)2)~2-Kr(sinh(/)2)^''], 
where we have varied the coefficient v over several orders 
of magnitude and taken r = 0, 1 in each case. Each po- 
tential is represented in Fig. 2 by a different curve, with 
the r = examples corresponding to the curves bound- 
ing the shaded area. We will not give further details here 
because, as we will show, they are not important. 

We find that the results for /at^ fall into two distinct 
regimes depending mainly on the duration of the conver- 
sion process. We measure the duration by specifying the 
number of e-folds ('Hubble times') by which the scale fac- 
tor a shrinks during the process of conversion; it is useful 
to distinguish the case of "rapid conversion," where the 
conversion lasts at most about 0.2 Hubble times, from 
the case of "gradual conversion," where the conversion 
lasts on the order of 1 Hubble time (The details of the 
reflection affect the direction of the trajectory afterwards, 
which has to be taken into account when incorporating 
this conversion mechanism within the full cosmological 
model; but this is beyond the scope of this paper.) 
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FIG. 2: The predicted value of the non-gaussianity parame- 
ter /jvz, as a function of the duration of the reflection. The 
shaded range represents results for the simplest ekpyrotic and 
reflection potentials allowing kz to range from -1-5 (top) to —5 
(bottom). The upper and lower boundaries can be extended 
(dashed and dotted curves) by carefully tuning the reflection 
potentials, but the change is not signiflcant for typical dura- 
tions 0(1). 



In the case of rapid conversion, the conversion is inefh- 
cient so that the leading order curvature contribution IZ^ 
is smaller by a factor of 0(10) or more compared to the 
case of gradual conversions. Reducing TZj^ by a factor of 
ten automatically amplifies all the contributions to /atl 
by a factor of 0(100), although the summing up of all 
the contributing terms cancels out some of the amplifica- 
tion. The net result for jjssL turns out to be of O(IOO) or 
more, as shown in Fig. 2 for the case when the duration 
of the conversion is only 0.2 Hubble times. The figure 
also shows that the result is sensitive to the form of the 
potential and the potential parameters (v and r) used to 
model the reflection. However, much of the entire range 
lies outside the observationally acceptable range, and so 
can be dropped from further consideration. Here we are 
chiefly interested in making general predictions, and not 
in the accidental cancellation of large numbers, which 
would clearly be flne-tuned and unattractive. 

The second regime is where the reflection and conver- 
sion process are gradual. In these cases, the conversion is 
efficient and all contributions to /atl tend to be reduced. 
Roughly, for K3 ranging from —5 to 5 and for typical val- 
ues of e ~ 50, we find that the predictions converge to a 
range 



'60 < / 



NL 



< 



-80, 



(19) 



with very weak dependence on the parameters of the re- 
fiection potential. If e is pushed to higher values, e.g. up 
to e = 200, this range expands to —120 < Jnl ^ -t-160. 
Analytic estimates of the above calculations confirm that 
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the results are rather insensitive to the specific form of 
the potential causing the reflection, but depend most 
strongly on the duration of conversion as well as roughly 
linearly on Kz^/e [24]. 

Several comments are in order. First, our results show 
that, in the ekpyrotic model, it is very unnatural to get 
a value of |/jvz,| ^ 1, the value obtained for slow-roll 
inflation [22]. The contributions to the local /jvl are 
generally greater than ten so that small values can only 
be achieved by accidental cancellation of large terms. 

Conversely, obtaining large local fNL is completely 
non-generic for inflation, requiring the finely-tuned ad- 
dition of extra fields and interactions. When the fields 
and interactions are introduced, there is no predictive 
convergence to a bounded range for /jvl- The result 
can be anything from negligible Jn l to values of either 
sign and extending orders of magnitude beyond current 
bounds. By contrast, in the ekpyrotic model, conver- 
sion of scalar fluctuations to curvature perturbations is 
absolutely required, and, at least for one proven con- 
version mechanism, generically generates a measurable 



fNL- We have found that a rapid conversion, whether in 
an ekpyrotic or kinetic energy dominated phase, gener- 
ically leads to a large I/jvlI — O(IOO) or greater that 
is already ruled out observationally. Furthermore, if the 
conversion occurs in the ekyprotic phase, it is difficult to 
obtain simultaneously an observationally acceptable /jvz, 
and an observationally acceptable spectral tilt. On the 
other hand, a wide range of ekpyrotic/cyclic models sat- 
isfying current constraints on the spectral tilt generate a 
measurable non-gaussianity consistent with current ob- 
servations if the conversion is gradual and occurs after 
the ekpyrotic phase is completed. 
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